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Abstract 

We reconsider a string theoretic inflationary model, where inflation is driven by n muhiple 
coincident Z)3-branes in the finite n hmit. We show that the finite n action can be continued 
to the hmit of large n, where it converges to the action for a wrapped D5-brane with n 
units of UiX) fiux. This provides an important consistency check of the scenario and allows 
for more control over certain back-reaction effects. We determine the most general form of 
the action for a specific sub-class of models and examine the non-relativistic limits of the 
theory where the branes move at speeds much less than the speed of light. The non-Abelian 
nature of the world-volume theory implies that the inflaton field is matrix valued and this 
results in modifications to the slow-roll parameters and Hubble-flow equations. A specific 
small field model of inflation is investigated where the branes move out of an AdS throat, 
and observational constraints are employed to place bounds on the background fluxes. 
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1 Introduction 



Modulo several signatures from the strong-coupling regime of QCD PQEIE]! as well as indirect 
evidence possibly emerging from terrestrial accelerators [57j, string theory remains best-tested 
in the realm of early universe cosmology. Observations of the cosmic microwave background 
(CMB) provide strong evidence for the fiduciary ACDM model of cosmology, together with 
a period of primordial inflation which set the initial conditions for the density fluctuations. 
Ideally primordial inflation occurred in an energy regime where degrees of freedom unique to 
string theory are excited; or, at the very least, string theory should provide the mechanism 
for an epoch of inflation that is consistent with current observations. An understanding of 
early-universe dynamics therefore provides a direct means of exploring aspects of string theory. 

Given this motivation, inflationary model building within string theory has become an estab- 
lished research direction, yielding many different models. Despite the large number of proposals, 
this pursuit serves — at the very least — to constrain the parameter space and overall viability of 
string theory [58]. A simple dichotomy of these scenarios emerges from the origin of the inflaton, 
i.e., whether this field originates from the open or closed sector of the theory. In both cases, 
numerous novel phenomena can arise and, in this sense, string theory provides a natural expla- 
nation for possible effects which can not be understood within the context of vanilla, slow-roll 
inflation. Since such results reside mostly within the framework of effective field theory, one 
cannot be certain if string theory is the unique UV completion. Nonetheless, it is remarkable 
that the classes of models selected by string theory have such a wealth of unique and potentially 
detectable phenomena. 

In this paper we focus on a popular form of brane inflation embedded in the open-string sector 
of string theory. This model of inflation associates the inflaton with a modulus parametrizing 
the separation of a probe brane from a stack of branes residing at the bottom of a warped 
throat. There has been immense excitement that such a model could lead to distinct observa- 
tional signatures, particularly regarding a non-vanishing bi-spectrum in the primordial curvature 
perturbation (see [47^ [48l [49l [50] for examples). However, there remain many open questions 
regarding the regime of validity of this scenario, at both the theoretical and observational levels 

[ElETllSl]. 

In view of this, we go beyond the simplest approximation and consider a more complicated 
model comprised of n multiple coincident branes that follow a single trajectory. Due to the 
nature of the world- volume theory living on these branes, the resulting field-theoretic analysis is 
automatically a non-Abelian theory, where the infiaton is now a scalar field transforming in an 
adjoint representation. One strong objection to using a theory of this form is that we wish the 
standard model to exist in another (warped) throat, and therefore a world-volume description 
of inflation is not desirable. However, we will argue that such a theory has a dual description 
in terms of a wrapped higher-dimensional brane, and therefore represents an important theory 
within the class of DBI-inflation models. Importantly, we will argue that the wrapped brane 
description is valid for large values of C/(l) flux by showing how the multi-brane action converges 
as we take the large n limit of the flnite n theory. 
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Since the fully relativistic theory has been (relatively) well explored [391 HOI HIl El, we 
focus in the present work on the non-relativistic limit, where the velocity of the inflaton field is 
significantly smaller than the speed of light. It is in this sense that we refer to such a regime as 
'non-relativistic' and this is distinct from the limit where the effective sound speed of inflaton 
field fluctuations is much less than that of light. In the standard DBI-inflation model, the 
non-relativistic limit of the theory is simply a canonical scalar field in the appropriate FRW- 
background. However, due to the non-commutative structure inherited from the non-Abelian 
field theory, the non-relativistic limit in our scenario is a non-canonical scalar field theory. 
Consequently, one should anticipate different physical effects when compared to the standard 
'slow roll' inflationary scenario. 

The paper is structured as follows. In section [21 we remind the reader of the string theoretic 
construction for this class of models. Firstly we argue that the large n limit of the multi- 
brane theory is precisely dual to a single wrapped DS-brane with n-units of U{1) flux. We 
then consider the finite n limit of the multi-brane theory, since this is analytically tractable, 
and proceed to argue that the finite n limit does indeed analytically converge with the large n 
theory in the correct scaling limit. This is the first time that such a result has been established 
in the literature. In section [31 we develop the framework for analyzing inflationary dynamics in 
the finite n theory, and construct a simple example of small-field inflation in an AdS background 
driven by an inverted harmonic oscillator potential. In section[4l we specialize to the n = 2 theory 
and consider in some detail how the inflationary observables differ from those of the single brane 
model. We then discuss how the flow equations for such a model can be constructed. Finally, 
we conclude in section [5] with some general remarks and a discussion of future directions. 



2 The Action for matrix cosmology 

2.1 Matrix cosmology at large n 

Numerous extensions of the vanilla DBI-inflation model have been proposed, yielding novel 
results and phenomena. These include, but are not limited to, multiple-field scenarios |19i [20l 
[2T1 [26] , multiple branes [221 ISHl [M| , monodromies [25] , Wilson lines [13 [l6] and the inclusion of 
bulk form fields [18]. In this section, we are interested in the more conservative generalization 
of the model, where the solitary Z)3-brane is replaced by a solitary D5-brane wrapping a non- 
trivial two-cycle of the internal geometry (usually taken to be an 5^ — note that this is a simply- 
connected space and therefore we must turn on some world-volume flux to stabilize the brane) 
[27\ [28] . Whilst there has been significant progress in the field of flux compactification in the 
type II theory onto (conformal) Calabi-Yau manifolds, we will consider a simpler class of warped 
background geometries based upon the well understood conifold geometry 

ds^ = h^ip)dsl + h-^{p){dp^ + p^dn]^^). (2.1) 

The basic idea is that fluxes in the compactification back-react to form a throat (parametrized 
by a radial coordinate p) over the base manifold, X^. We allow the brane to be localized (and 
flat) in the large dimensions, wrapping a two-cycle within the internal manifold. We further 
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assume that the brane is dynamical along the throat direction and freeze out any angular degrees 
of freedom. 



Without specializing to a particular supergravity background solution for the warp factor, 
h{p), we know that generically one must turn on magnetic flux along the two-cycle directions in 
order to stabilize the configuration. The U{1) field strength must be proportional to the volume 
of the wrapped cycle, where the constant of proportionality is determined by the n units of fiux 
threaded through the 5^: 

F^'^ = ^^2, (2.2) 
where uj2 is the associated two-form of the sphere. 



The corresponding action for a DS-brane is then given by the usual DBI expression (see [29j 
for further discussion): 

S = -nj d6xe-*^-det(Gab + B^b + AF^b) - J C^e^ j e^^, (2.3) 

where a "denotes the pullback of that particular tensor to the world- volume, and A is the inverse 
of the fundamental string tension which couples to the U{1) field strength. The field ^> is the 
dilaton defining the string coupling constant in the low energy theory. However, since we are 
assuming that the throat is sourced by D3-brane charge, we can set this term to unity without 
loss of generality. For more involved backgrounds, the dilaton will indeed be non-trivial and can 
lead to more complicated dynamics. 



The second term in Eq. (|2.3p corresponds to the coupling of the brane to the bulk form 
fields in the RR-sector once they are pulled-back to the world-volume. The presence of the 
summation indicates that there can be coupling to form fields of lower degree, provided that 
there is a non-zero B (or F) term. Note that T5 and /is denote, respectively, the tension and 
charge carried by the brane, and that these are related via supersymmetry. By utilizing the 
above metric and gauge field ansatz we can compute the action for the wrapped Z)5-brane. It 
is given by 

S = -47rT5 J (fi ^/iVl - h-^p^^jhi^X^n'^ + / - /i^^A^ , (2.4) 

where we have worked in physical coordinates and integrated out the directions along the S^. 
For cosmologically relevant solutions one must minimally couple this term to the usual Einstein- 
Hilbert action. For a concrete embedding of this action into a particular string theory back- 
ground we refer the reader to [28]. 



One of the many interesting features of string theory is the existence of duality symmetries, 
whereby a brane configuration can be related to another (equivalent) one. In this instance, we 
can employ our knowledge of the Myers dielectric effect [29| to understand how the above action 
is captured in the dual picture by n coincident Z)3-branes expanding along a fuzzy two-spher^. 
The full discussion is intricate and chronicled at length elsewhere (see \2>Q\ i34j for examples in 
non-trivial backgrounds). We therefore quote the results and refer the reader to [40^ I41j for 

''See the alternate proposal by [ST] . 
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more details on the relation to cosmological model building. The validity of the Myers action 
using boundary fermions has been considered at length in [321 [33] , and can certainly be trusted 
at leading order. An important related usage of the action is discussed in [35]. 

The action for n coincident branes in the limit where n ^ 1 can be written as 

S = -nT, J d^h^ (^i-h~^p2^i + ^^^- 1^ , (2.5) 

and it can be shown that this is exactly the same as the action (j2.4p abov^ by employing the 
known relationship between the brane tensions, = An'^a'T^, and by identifying the large n 
limit in both cases. This implies that in the (macroscopic) D5-brane theory, we must take the 
flux to be large. Since our background is relatively simple, one could also consider the S-dual 
theory, where the Z)3-branes expand into a wrapped NS5-brane. In this case, the question of 
where the standard model degrees of freedom reside becomes crucial since fundamental strings 
cannot end on the world-volume of the fivebrane. Within this context, a fivebrane model was 
recently proposed as a concrete realization of axionic monodromy inflation [59]. 

An important result is that the back-reaction of higher-dimensional branes becomes impor- 
tant in the relativistic regime [27]. Given our dual interpretation of the L>5-brane action, this is 
intuitively obvious, since we are essentially taking n ^ 1, and this will typically back-react on 
the geometry and invalidate the probe-brane approximation. Thus, it is clearly desirable to con- 
struct a theory where n is not taken to be extremely large. This should reduce the back-reaction 
and allow for more control over the theory. However, constructing the action for multiple D3- 
branes in the finite n limit is a highly technical issue and has only been resolved in a few simple 
cases. Nonetheless, progress has been made and in the following subsection, the key results that 
are relevant to the present discussion are summarized. Full details of the construction can be 
found in 1391. 



2.2 Matrix cosmology at finite n 

The tractability of the problem for finite n is intimately tied to the symmetrized trace prescrip- 
tion (STr) associated with open string scattering amplitudes. This suggests that one must first 
average over all symmetric permutations before taking the gauge trace. Fortunately, one can 
make headway by identifying the scalar fields with a finite dimensional irrep of S0{3) ~ SU{2), 
which leads to the following result [361 EZ] • 

n/2 

STr{a'ay = 2{2q + 1) ^(2i - l)^? n even 

i=l 

(n-l)/2 

= 2(2g + l) ^ (2i)2'? n odd. (2.6) 

i=l 

^At least to leading order. Examining the duality to higher orders is certainly interesting, but highly non-trivial. 
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As an aside we remark that for the flat D3-branes, the scalar fields parametrize the direction 
of a transverse S*^ — assuming that the full, non-compact, background geometry admits a sub- 
manifold with an S0{3) isometry. The corresponding action for n multiple L'S-branes is then 
given by the usual expansion of the Myers DBI action: 

P„ = -T3STY lh\p) f; {-ZR^)''YP(a'a')''+p(^^j^^^(^^^^^+V{p)ln-h\p)lA (2.7) 

\ k,p=0 j 

where we have defined the following terms 

Z = X'h-Hp), Y = AX'R^h-Hp), (^^^) = r (2.8) 



qj T{3/2-q)T{l + q) 

and the fuzzy sphere radius R is related to the physical radius p and the number of branes via 

p^ = X^R\n-lf. (2.9) 

Those familiar with such a construction should note that one usually relates the physical coor- 
dinate to the quadratic Casimir (n^ — 1). However, for finite n there is a more precise definition 
of the physical radius in terms of ratios of operators [37]. This is a necessary requirement for 
the convergence of the theory, and agrees with the definition in the limit of large n. It should be 
noted that the radial coordinate p, the potential V(p), and the warp factor h{p) are all singlets 
under the symmetrized trace prescription. Moreover, the potential as written here is dimension- 
less, since we have absorbed a factor of the brane tension into its definition. It is also important 
to note that the physical radius is a function of n in this model and this will play a role when 
the recursion relations are employed. The key point is that the scalar field, which we wish to 
identify with the inflaton, is now matrix valueqj and is in a finite dimensional representation of 
the corresponding gauge group. 

The most important representation of SU{2) is the fundamental (or two dimensional) one, 
since the representation is simply that of the Pauli spin matrices. It also turns out that the 
form of the symmetrized trace ansatz leads to a recursive structure for the theory, where we can 
construct the action for any n once the action for n = 2 has been computed. Thus, the correct 
way to write the action is in terms of spin-^ variables as first elucidated in [39]. It is found that 
the action takes the form 



(n-l+<5„)/2 

PniZ, ¥)= Yl ^2((2A; - SnfZ, {2k - 5nfY) - nT^{V{p) - h\p)) 



(2.10) 



where <5„ = 1 when n is even and 5„ = when n is odd. A similar recursive structure exists for 
the energy density, En{Z,Y) (more specifically the time-time component of the kinetic part of 
the energy- momentum tensor). 



^For recent additional work using the same idea we refer the interested reader to [56] . 
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For completeness, we write down the recursion functions P2 and E2- 
P2{Z,Y) 

E2{Z,Y) 

where Z, Y are the functions defined in Eq. (j2.8p . Note that the algebraic forms of these 
functions are significantly different from those of the n = 1 action and, consequently, one should 
anticipate different physics to emerge. 



2T3h^ jl + 2Y -{2 + 3Y)ZR'^ 
~VT+Y \ Vl - ZR^ 

2Tsh'^ fl + 2Y -YZR^\ 

VT+Y \ (l-Zi?2)3/2 1 



2.3 Convergence with the large n hmit 



A crucial question that now arises is whether the finite n action defined in Eq. ()2.10p really 
does reproduce the action (j2.5p in the large n limit. We present here, for the first time, an 
argument which suggests that this is indeed the case. Our approach is to expand both actions 
as Taylor series and compare coefficients term by term. We focus on the NS-NS sector, since 
the -R-R-sector converges in a trivial manner. To proceed, we rewrite action ()2.5p at large n in 
the more compact form 



P 



1 + ^ 

C 



1/2 



where we have introduced the following variables: 



1 - (n- 1)2^/^2 



(2.12) 



In 



Trig 

\h 



C = n^ 



Y 



In 



{n-iy 



(2.13) 



The canonical inflaton field is defined in terms of the throat geometry via the standard relation 

4> = Vt^p- 



Expanding the DBI action requires the introduction of the binomial coefficient, which we 
write as 



(l + x) 



00 

E 

j=0 



a 



a 



a{a — 1) . . . (a — j + 1) 



(2.14) 



For the remainder of this subsection we explicitly drop the summation sign, although the sum- 
mation over Latin indices is always implied. The key point is to observe that the ratio Yq/C 
will generally be small in the large n limit. This implies that we can expand the action ()2.12p 
up to the leading order terms in {ip'^yY'^. It follows, upon expansion, that 




(2.15) 
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The action (j2.15p is what we aim to reconstruct using the finite n formahsm. We demonstrate 
this exphcitly for odd n (the case of even n is analogous). The relevant n = 2 term, written in 
terms of the new variables (j2.13p . is given by 



2r3/i4 



i + 2y-(2 + 3y)V'^) [i-Tp 



-1/2 



(2.16) 



We now proceed to expand the velocity factor (1 — tp'^) ^1"^ up to terms of order (ip'^y ■ After 
some algebra, we find that the relevant term in the expansion is given by 



P2 ^ 



VT+Y 



1 + 2y - (2 + sy)^'^ 



+ 



(i - 1)! 



(j - 1)! 



(1 + 2j) ^ (2 + 2j)Y 



2j 



2j 



(2.17) 



Performing the analogous expansion on the (1 + Y) factor up to terms of order then 
yields the result: 



P-? 



2r3/i^(l + 2j + 2i) r/ 1\ f3-2i 
4i!j! [V 2 

-2nh^ f l±^f±^') QyQ^ .... 



3-2j 



i-r) 



2V 



(2.18) 



However, we should recall that this is simply the expansion of the n = 2 action and we must 
employ the recursion relations (j2.10p in order to generate the full structure for any value of n. 
For odd n, this is achieved through the rescaling 



tP'^ {2kfi)'\ Y {2kfY 



(2.19) 



Furthermore, we must also introduce a new summation over the k variable to obtain the complete 
expression for the leading order terms in {ip'^yY^ in the action P„. Taking these considerations 
into account, we arrive at the expression 



1 + 2j + 2i 



(n-l)/2 



(2.20) 



k=l 



To establish the correspondence with the large n action (j2.15p . we need to trade the sum over 
k for a function of n. This is achieved by considering the well-known algebraic condition: 



n' 



P+i 



p+1 



< 



k=l 



< 



{n + l)P+i 
p+l 



(2.21) 



It can be seen that in the limit of large n the lower and upper bounds converge. Thus, we can 
make the following algebraic identification 



V- n2(i+j)j^2(t+j) _ [n- L)^^ 

^ 2(l + 2i + 2i) 



(2.22) 
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Finally, inserting this relation into the action (j2.20p . collecting together all the terms and sub- 
stituting for Y in terms of Yq, leads us to the result 




Eq. (|2.23|) may be compared directly to the action (|2.15|) in the large n limit. 



(2.23) 



In conclusion, therefore, we have argued that the action for finite n does indeed converge to 
the full action at large n, which is a non-trivial result. This indicates that the finite n action can 
be trusted (at least for the simple case we consider) . Whilst the full structure of the non-Abelian 
DBI action remains unknown, the Myers action [29] — and therefore our finite n theory given by 
Eq. (j2.10p — will be a good approximation to leading order in an a' expansion. In the following 
section, we consider some of the consequences of this action for the infiationary scenario of the 
early universe. 



3 Matrix inflation at flnite n 



3.1 General remarks 



One of the main characteristics of DBI-driven inflation is that the infiaton can move relativis- 
tically and still drive a sustained period of accelerated expansion. This is possible because the 
warped metric redshifts all physical scales associated with the brane dynamics. This feature is 
frequently exploited to simplify the functional form of the solution by considering the so-called 
'ultra-relativistic limit', where the kinetic contribution to the action is rewritten in terms of a 
function, 7, corresponding to a generalization of the relativistic factor of special relativity. 

In our scenario, the algebraic structure of the multi-brane action forces us to consider a more 
general version of this function, which we parametrize in terms of the inflaton (j) such that 

7.= (1-T^(^i-^) 1 ^ (3.1) 




Note that there is an implicit summation over k in this expression. Let us study this function 
on its own for the moment. It is trivial to see that increasing the number of branes, n, increases 
the number of poles of the function. For example, we see that in the case of n = 7 



-1/2 / X -1/2 / .„ X -1/2 



''=y-ih) ^['-^ih) 

and it follows that 7 diverges in the same limit as that of the usual DBI-inflation models, namely 
when (])'^ ~ h'^T^. The divergences arising from the higher order expansion in k are not physical. 
Consequently, it is the lowest term in k that determines the sound speed for infiaton field 
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fluctuations, at least if the analysis is limited to the behavior of the 7-function. The implication 
is that 7fc ~ 7 in the ultra-relativistic regime. 



If we include the coupling to Einstein-Hilbert gravity, we find that the energy density of the 
DBI-scalar is not conserved, due to the presence of the FRW scale factor. It is a standard result 
that En = —3H{Pn +En), where H is the usual Hubble factor, and from this expression we can 
then write down the dynamics of the Hubble parameter as follows: 

Hn = -m{P + E) (3.3) 

This is a highly complicated sum over all k states. However, the analysis can be simplified 
somewhat if it is assumed that the generalized 7-function p.ip is large, i.e. 7^ ^> 1 V A:. If this 
condition is satisfied, we find that 

P k 

where we have explicitly included the dependence on the number of branes in the H term 
and used the fact that 7^ ~ 7 in this regime. This expression differs significantly from the 
corresponding limit of the single brane (n = 1) solution. 



3.2 Non-relativistic expansion 



What happens to the expansion of the action in the non-relativistic limit? Related work for the 
single brane case is discussed in [541 [55] . Fortunately, one can see that the symmetrized trace 
operation commutes with the non-relativistic limit and therefore we can immediately write down 
the following equation of motion for the inflaton: 





Here a prime denotes differentiation with respect to (p and we have assumed the validity of the 
Hamilton-Jacobi expansion when higher order terms in the velocity expansion are neglected. 
Given an expression for the time dependence of the inflaton, the corresponding slow-roll param- 
eter e = H /H"^ can also be deduced: 




Vl + Y{2k - 5nf 



(3.6) 



It follows, therefore, that the slow-roll parameter is suppressed relative to that of canonical 
inflation by the additional sum in the denominator, which we will denote by for simplicity, 

i.e.. 



Ml /H' 



p 



- 



(3.7) 
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The maximal suppression occurs when J2k fk takes its largest value, although one must also be 
aware that Y is itself a function of n, which complicates the analysis. Indeed, it is trivial to see 
that 

and so the dependence on the inflaton is essential in determining the dynamics. 



In most cases of interest, the warped throats fall into one of two classes. The first is an 
AdS-iype throat, arising for example from the near horizon limit of D3-branes. In this instance 
the warp factor takes the form /i ~ </> and Y is fixed at some non-negative value. The second 
class corresponds to regularized throats, where the back-reaction of the form fluxes are used to 
'cap' the throat at the tip. A canonical example of such a background is given by the warped 
deformed conifold (BJ [7j. The warping in this case tends to a constant at some finite value of 
the throat, and therefore we find that Y is an increasing (decreasing) function of the inflaton 
depending upon whether we are considering small (large) field inflation. For the IR scenario 
(an example of small field inflation), one sees that fk is rapidly dominated by smaller values 
of n and therefore inflationary trajectories for n = 2 are preferred over those with larger n. In 
the UV (large field) case, on the other hand, we find that this function is sensitive to the precise 
value of the constants. Consequently, we cannot identify a priori whether more suppression 
occurs for smaller values of n, although one should remember that UV inflation is more severely 
constrained than the IR scenario. 



In general, we can write the Hubble parameter in the following form: 

3H^M^ ~ 2Ts Y,Wk{l + Sk) (3.9) 



k 

where the generalized effective potential is given by 



Wk ~ h\l + 2Y{2k - 5nf) + ( "^-T-^ {V - h^Ul + Y{2k-5n? 

\n-l + dnj 



(3.10) 



and the parameter 



{3 + 4Y{2k-5^f) 

~ 2T,{n - lyWk ^l + Y{2k-5nY 

quantifies the ratio of the kinetic and potential energies. It follows, therefore, that the generalized 
master (Friedmann) equation for the Hubble parameter can be expressed in the form 

2 



V- (3 + AY{2k - 5n? H'^M^ ( I \ , , 

^ ^l + Y{2k-5n? {n-l? [EpfpJ ^ ' 

and this equation can be employed to determine the inflationary trajectories. 



3.3 A matrix inflation model in AdS 



In this subsection we investigate a concrete example of non-relativistic matrix inflation with the 
aim of identifying regions of parameter space that are consistent with cosmological observations. 
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Specifically, we consider IR (small field) inflation driven by an inverted harmonic oscillator 
potential, where the branes are propagating out of an AdS throat. 



It is simple to show by expanding Eq. (j2.10p that the leading order expression for the 
pressure can be written as 

P,^Q,^-nT,Vi,Hn.'(n-E^02l^] ,3.13) 



where we have defined 



^^ (3 + 4(2fc - 6nrY){2k - 6r,y 
Y (n - 1) Vl + - ^„)2y ' 



For a class of AdS backgrounds, the warp factor h{(j)) = 0/ {L^/T^) is determined in terms of the 
AdS radius of curvature L = [Air^ggM K / (Vol{X^)mg)]^^^ , where gs denotes the string coupling, 
Vol(X5) is the volume of the base and the background fluxes M,K thread through the three- 
cycles in the compactification. Eq. (jS.Sp then implies that Q is independent of 0. We may 
therefore introduce a new field, ■0 = V^Q4'i order to write the pressure in a canonical form. 
The effective potential energy of this field is then given by 

2iTMKrr-\ ^ ^1 + (2k - S^r-Y J 

where we have written the explicit parameters 

Since typically the product of fluxes must satisfy MK ^ 1, it follows that the second term in the 
effective potential is suppressed and therefore only the term proportional to V{ip) is expected 
to dominate the dynamics. 

Indeed, if one assumes that the inflaton potential takes the form of a simple inverted harmonic 
oscillator, V{(j)) = nT^VQ{\ — I54P' /M^), where /? is an order one parameter of the theory, we may 
write ^ ^ 

W{i,) = 1^0 - ^ (3.17) 

where 

Wo = nT,V,, -^ = ^- (3.18) 

Coincident branes located near the tip of a warped throat are expected to be attracted towards 
branes residing in other throats through the generation of such a tachyonic potential. Since 
our theory is now in canonical form, we can adapt the usual inflationary tools to our current 
objective. During inflation the Gaussian curvature perturbation has an amplitude given by 

and spectral index 

1 - = 6e - 2Ty (3.20) 
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where e = ^Mp{W' /W)"^ and rj = M^W" /W are the usual slow roll parameters obtained from 
derivatives of the effective potential. Near the maximum of the potential we see that 

which implies that e ^ jr^l ^ 1. The spectral index of the density perturbations can therefore 
be computed as 

1 - n, = ^ (3.22) 

This has an interesting dependence on the number of coincident branes, since Q is not monotonic. 
Indeed, Q has a local maximum when n takes the lowest possible value (n = 2 for the even case), 
which suggests that 1 — is initially small. As we increase the number of branes, we see that 
1 — 71^ also increases until it attains a local maximum around n ~ 6 (for the even case), before 
asymptotically tending to zero as n increases further. 

The above expressions assume implicitly that the slow-roll conditions hold until the very 
end of inflation, i.e., it is assumed that the branes continue to move non-relativistically. The 
self-consistency of this assumption can be verified by employing the effective field equation 
?,Htp ~ -dW/d^ and the Friedmann equation SiJ^ ~ Wq/M'^, together with Eq. The 
non-relativistic limit, cjp' <C T^h^, is then equivalent to 

^2»/3L^r3(l-n,)^ (3.23) 

Since the inflaton is a monotonically increasing function in this scenario, it suffices to show that 
the bound (|3.23p was satisfied when observable scales crossed the Hubble radius. The value of 
the field at that time is related to the normalization of the CMB power spectrum such that 

V'cmb _ 1 nT,Vo 

and substituting Eq. (|3.24p into the constraint (|3.23|) then yields the consistency condition 

pTsL^l - nsfn^Vn « 1 (3.25) 
This constraint can typically be satisfied for a large region of the physical parameter space. 

The value of V'cmb can also be determined from the scalar field equation of motion to be 

V'cmb = V'ende-^^-"^)^/' (3.26) 

This allows us to write the tensor-scalar ratio in the following form 

r = 16e = 2(l-n,)2^ (3.27) 

p 

suggesting that gravitational waves are negligible in this model. 

We may gain further insight by specializing to the case of even or odd n. Let us choose the 
latter case for clarity and reconsider the constraint on the spectral index (|3.22|) . We can clearly 
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see that since Y is positive-definite, the definition of Q given in Eq. ()3.14p imphes that this 
parameter is bounded from above such that 

Utihzing the bound (j2.21|) enables us to write down a further constraint that is shghtly weaker, 
but expressible simply as a function of n: 

where we neglect the terms coming from higher orders in the expansion. 

The WMAP five-year data implies that Ug > 0.93 {2a) when the gravitational wave back- 
ground is negligible. This in turn imposes the condition Q > 30/3 ^ n if we take (3 ~ 
Now, since the first term in the expansion of (j3.29p is essentially 0{n), and therefore negligible 
compared to Q, the leading order bound is dominated by 

64^0 fn + l\^ , , 

and this results in Yq being bounded from below: 

y 5/3 (n - If 

- (l-n.)(n+l)5 ^'-''^ 

In principle, therefore, condition ()3.3ip may be interpreted as an observational lower bound on 
the product of the fluxes MK after 'canonical' values for the string coupling, Qs — 10~^, and 
base manifold volume, Vol(X5) ~ tt^, have been specified. It follows from Eq. (13.160 that 

- (1-n,) (n + l)5 ^ ' 

Moreover, the ratio (n — 1)^/ (n + 1)^ is a monotonically increasing function of n and this implies 
that consistency with observations would be difficult to achieve if condition ()3.32p was not 
satisfied for n = 3. For example, invoking the central value for the spectral index inferred from 
WMAP, Us = 0.96, and setting (i ~ 0(1) suggests that MK > 630. One should note that the 
fluxes depend on the inverses of both the string coupling and the spectral index, and a more 
weakly coupled theory therefore leads to a tighter bound on the fluxes. Likewise, it becomes 
more difficult to satisfy this condition as the spectrum becomes closer to scale invariance. 

To summarize, in this model the spectral index of the density perturbations is determined 
by the parameter Q that quantifies the non-Abelian nature of the multi-brane configuration. 
Consistency with observations is possible if the background fiuxes are sufficiently large. 

In the following section, we proceed to investigate the dynamics of the n = 2 configuration 
in more detail. 
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4 Cosmology of the two-dimensional representation 



Given the recursive structure of the action in Eq. (|2.1U|) . we see that the two dimensional 
representatioE0 of SU{2) is the most important, corresponding in our matrix language to a 
configuration of two coincident L)3-branes. Note that this is not the same model as that con- 
sidered in [23], which considered two D3-branes separated by a distance larger than the string 
length. The world- volume gauge theory in that case is therefore U{1) x U{1), whereas our 
model possesses a U{2) symmetry. Note that this is the gauge group of the open string states on 
the world-volume, and should not be confused with the SO{3) ~ SU{2) symmetry group that 
parametrizes the transverse space. 



4.1 The relativistic limit 



To begin this section let us focus on the relativistic limit to highlight how the theory differs from 
the n = 1 solution. Unlike the case of a single Z)3-brane, the pressure and energy of the model 
have very different dependence upon the brane velocity. Given our generalized definition of 7^, 
it is convenient to now introduce the following notation 

72= , ^ ■ ■ (4.1) 
VI - ZR^ 

As in the case of n = 1 the sound speed of the Fourier modes, c^, is reduced from unity — the 
relevant calculation proceeds in the usual manner and we find that 



1 / 27|(3 + 4y)-(7|-l)(4 + y) , 

27iV 7l(3 + 4y)-(7i-i)y ' ^^-'^ 



which is very different from that of the standard DBI inflation models (see [H] for example) . In 
the ultra-relativistic limit the above expression simplifies to become 



cl^^A-f:^] (4-3) 



2 _ J_ ( 2 + 7Y 
" ^ 67I \ l + Y 

and this should be compared to the ~ 1/7^ dependence of the n = 1 model. One sees that for 
y S> 1 the sound speed varies as (7/6)7^^, whilst in the converse limit it behaves as (1/3)7^^. 
More importantly for cosmology, the sound speed modifies the tensor-scalar ratio as follows: 

(4.4) 

with the presence of an additional enhancement factor due to the contribution of Y. Indeed, 
this accounts for the enhancement of r, relative to single brane models, discussed in |39j . 

Given the modified expression for the sound speed, and the non-linear form of the action, we 
can compute the bi-spectrum arising from higher order interactions of the different Fourier modes 




^Note that this is the adjoint representation of the group. 
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in the density perturbations. For simplicity, we focus on the equilateral triangle configuration to 
highlight the relative contribution to the non-Gaussian parameter fni- A standard calculation 
reveals that 



_857| / 7i(3 + 4y)-y(7|-l) \ 8^ 
162 \2^l(3 + 4Y)-i^l-l){4 + Y)J 324 

5(7|-i) /7l(5 + 6y)-y(7|-i)' 



81 V7|(3 + 4y)-y(7|-2) 



(4.5) 



This function is negative definite for all values of 72 and Y in the relativistic limit. Moreover, one 
sees that the detectable range of fni falls into the regime where y <C 1. For example, in AdS-type 
backgrounds, where h ^ (j), the parameter Y reduces to a constant determined by Eq. ()3.16p . 
For sufficiently large MK/g^, we find that Y is suppressed, and detectable levels of (negative) 
fnl are therefore more likely. Note that in the limit where y — > 0, the functional form of the 
action becomes very similar to that of the single brane models. Thus, a small (but manifestly 
non-zero) value of Y is crucial for generating a detectable signature of non-Gaussianity in the 
primordial curvature perturbation. 

Although the highly non-linear nature of the multi-brane action leads to interesting observ- 
ables, it is also of interest to investigate the nature of the inflationary trajectories. This can be 
achieved by calculating the variation of the Hubble parameter from Eq. (|3.3p . which takes the 
form 

After taking the ultra-relativistic limit of Eq. (j4.6p . as well as the corresponding limit in the 
definition of the energy density 

E2 ~ 2r3 (hSlVl + Y + V - /i^) (4.7) 
we see that the 'slow roll' parameter, e = H / H'^ , becomes 

e^^fl + ^^l^V7l + lil±^^ (48) 

This leads to a novel phenomenon in multi-brane DBI models. The final term in parenthesis is 
essentially a term of order one suppressed by the relativistic factor and it effectively decouples 
in the relativistic limit. If the effective potential iV — /i^) is much smaller than the kinetic term, 
inflation will never begin. On the other hand, if the potential term dominates — even for the fast 
rolling solution — we see that this is a concrete realization of eternal inflation. 



4.2 The non-relativistic limit 



We now proceed to consider the non-relativistic limit of the theory. In this limit, the energy 
density takes the following, almost canonical, form: 

ii;2 = ^== (vF(c/>) + ^(3 + 4y)j (4.9) 
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where the effective potential is given by 

W{^) = h^{l + 2Y -VTTY) + V{(I))VTTY. (4.10) 
We also note that the pressure, P2, in the non-relativistic limit is given by 

P,{Z, Y) = -^== {w{^) - ^(3 + AY)^ . (4.11) 

As an aside, we briefly comment on the non-Gaussianities in this limit. In the n = 1 case, the 
non-relativistic limit reduces to a canonical scalar field theory, and therefore one expects that 
will be of the same order as the slow roll parameters. For n = 2, it is convenient to parametrize 
the non-relativistic expansion as 72 ~ 1 + ^ and focus on the correction term 

5^ /210 + 167y\ ,^ 
fni ~ [ o , ... J (4-12) 



234 V 3 -h AY 

One can see that this term is very small even in the limit where y 3> 1. 

Assuming the validity of the Hamilton-Jacobi approximation, we can use the continuity 
equation to obtain the following expression for the time-dependence of the scalar field 

which is modified from that of the canonical theory by the y-dependent terms. Interestingly, 
we can write the slow roll parameter e in in exactly the same way as the canonical model 

3e2 



1 + €2 



(4.14) 



so that 62 < 1/2 is required in order for inflation to occur. Recall that £2 is the ratio of the 
kinetic energy to the potential energy in the theory and follows from the generalized definition 
in (|3ll]) . 

The final step for our analysis is to calculate the master (Priedmann) equation for the Hubble 
parameter. This is given by 

which follows from the generalized expressions derived in section 3.2. Since an inflating solution 
requires that the terms on the right hand side are smaller than those on the left hand side, we 
can iteratively expand the Hubble parameter in a derivative series. The decoupled leading order 
solutions are given by 

~ ?,mIhI-2WT^ (4.16) 

(3 + 4y) 



^^^ + Y 

(3 + 4y) 



Ml{H[^ + 2H',H',) ~ m,H2^^^^ 
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for the first few terms in the expansion and this enables us to reconstruct the Hubble parameter 
at leading order as follows: 




^ 1+^1 ' .^ +--- ■ (4.17) 



The validity of this expression requires that the second order terms are negligible in comparison, 
which is equivalent to the condition 

^ ^/T+Y /W" W'Y' [5 + AY) MpW I 2 \ 

^ 2 {3 + 4Y)\W~W^~ 2W (l + y)(3 + 4y)~ 12W y ST^W J ' ^' ^ 

One can see that in the asymptotic regimes for Y, this constraint is essentially a bound on the 
values of W and its derivatives, exactly as in canonical slow roll inflation. This is particularly 
clear for the case where Y is constant, since this implies that the effective potential depends 
solely upon the derivatives of the warp factor h{(j)) and the scalar potential V{4')- The limit 
where y — > is essentially that of the usual DBI scenario (modulo a field redefinition) and 
we can therefore be sure that inflation will occur. In the limit of large Y, we see that the 
higher-order corrections to the Hubble parameter become arranged into a 1/1" expansion and 
are therefore suppressed relative to the leading order terms. 



4.3 Dynamical trajectories 



In this section, we are interested in studying the flow equations for the n = 2 solution along the 
research- lines of \A2\ H3| IH] in both the relativistic and non-relativistic limits. Recall that in 
the case of a single brane, the inflaton equation of motion (in the Hamilton- Jacobi formalism) 
is given by 



2M^H' 



7 

and the flrst slow roll parameter en by 



(4.19) 



In the case of flnite n, the non-linear (and occasionally higher derivative) nature of the action 
prevents us from analyzing the flow equations in the usual manner. However, there should exist 
regions of parameter space where one can work with equivalent algebraic structures. As such, 
we may extend the relations flrst established in [H] by introducing the following generalized 
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flow parameters: 



2M^ fH'^ ^ 



/(0) 
2Mg H" 

2M^H'f{ct>y 

fi4>) H fi<p) 



(4.21) 



'A 



'a 



2Ml\ 










I 







1 S+^H 



i-i 



1 

/(0) ^i'+V 



where f{4>) represents a generalized function that needs to be determined in a limit of the theory 
where higher derivative terms can indeed be neglectecif]. 

It is convenient for our discussion to note the following relations between the parameters: 

r] = ^X, C = ^A, p = ^a, a = '^a. (4.22) 

It then follows that the dynamics of these parameters can be established by solving the set of 
differential equations 



de 

drj 
dNe 

dK 
dNl 

dp 
dNl 
S\ 
dNe 
d^a 
dN. 



-e(2e - 2ry + k) 

-K{2K + e-T]) + ep (4.23) 
-2pK + a 

-^X{lK + k-ri{l-l))+^+^X 
-^a{{l + 1)k + l)e - (/ - 1)7]) + 



a, 



where iVg represents the number of e-foldings before the end of inflation and is defined by 
the condition dNe = —Hdt. The critical (stable) points of the phase space trajectories then 
correspond to the vanishing of each term on the right hand sides. 

A simple calculation shows that at leading order in the derivative expansion, the function / 
for the n = 2 case becomes 

/~^M=. (4.24) 



*The flow parameters used in this context are only applicable to the leading order dimension four terms. One 
can include the higher derivative operators, but the flow equations cannot then be written as simple functions of 
{H' /H). We hope to return to this more general problem in a future publication. 
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This is significantly different to tlie analogous term for the single brane model, where / ~ 7- 
Moreover, since Y can also depend on the field 4> (provided h is constant), we find that there is an 
interesting scenario where the branes are moving relativistically, but / ~ constant. This feature 
will be important in understanding the full dynamical phase space flow. In the non-relativistic 
limit we can truncate the theory such that 72 ~ 1, and then we find that / is only a function of 
Y. In this case, / can be constant if the warping is of the AdS form h ^ (p. Thus there are two 
solution branches that are of interest if we wish / to be constant: 

• Relativistic limit and taking h ~ constant. 

• Non-relativistic limit and taking h ^ (p. 



Note that the latter branch leads to the strongest constraint, since / is always constant if the 
warping is of the AdS-type. Indeed, in this case we immediately see that k = = 'a which 
reduces the dimension of the overall phase space. Since the flow parameter e is first order, we 
take this to be the fundamental parameter for our analysis and outline the results below: 



e = 

Setting e to be zero fixes H' to be zero and therefore the scale factor can be trivially 
integrated to yield 

a{t) ~ e^* (4.25) 

Here we are forced to choose solution slices where rj = e and also 'A = e', as can be seen 
from the fiow equations. Again from this we can infer that H is fixed to be exponential in 
form 

where we must take the positive sign above and also assume that (/> is a decreasing function 
if this is to lead to an inflationary trajectory. From this we can then reconstruct the solution 
for the infiaton field velocity with the result that Ht ~ 1/e. Therefore we obtain the usual 
power law trajectory with a ~ t^^'' as in the literature 



These are the only possible behaviors in this regime. The final case to consider is when / is no 
longer constant but scales in such a way such that e, k remain constant. In this instance we can 
solve most of the dynamical equations exactly, obtaining relations such as 

?? = i(2e + K). (4.27) 

The key point is that one can reconstruct / in terms of the flow parameters. It will be convenient 
to perform the analysis on a case by case basis. Let us first consider the non-relativistic regime 
of the theory, where / ~ + Y . We note immediately that /i ~ (/> yields constant Y and 

therefore constant /, forcing us to consider the situation where h ^ constant. We now have two 
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choices depending upon whether we wish to consider smaU or large field inflation. In the former 
case it can be seen that f ^ as the branes move and therefore we obtain 



i2 



fa ~8M^./. (4.28) 

where is a constant of integration. It can be verified that / ^ as ^ diverges. For large field 
inflation, we can write the function as follows: 

which satisfies the required boundary conditions. The Hubble parameter can then be recon- 
structed trivially (in both cases) to yield 



H^Hoe^pi^. —r y^dcP (4.30) 




Since the resulting expressions are rather complicated, we do not write them explicitly. We 
note, however, that there are many different choices for inflating trajectories depending on the 
particular choice of constants (and their signs). 

The relativistic case is more complicated, since / can still remain constant if the various 
parameters scale in the right way. For simplicity, let us consider the case where Y is constant 
and we are also in the large field branch of solution space. This requires us to consider a limit 
where 72 diverges when — > in such a way that we can re-arrange the equation to solve for 
the relativistic factor as a function of the inflaton 

This yields the result that cfP' oc exp(— S^iV) and implies that the field falls off significantly faster 
than in standard single brane inflation. After substituting this dependence into the definition 
of the flow equation for e, we find that 



K 



H^r^, 5 = i, Hii±22 (4.32) 



In this case, therefore, the inflationary dynamics are sensitive to the ratio e/n, as in the single 
brane case, but also to the precise value of Y which is, in turn, set by the value of the background 
throat charge. 



5 Discussion 



In this work wc have considered the cosmological consequences of using the action for n coincident 
-D3-branes in the finite n limit. We have argued that the action for n coincident -D3-branes, with 
n S> 1 and with scalars transforming under the n dimensional representation of SO (3), is the 
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same as that arising from a single D5-brane wrapping an 5^ and carrying n units of U{1) flux. 
Through the usual notion of string duality, this implies that we can either consider cosmology 
on n coincident D3-branes, or a single wrapped DS-brane with commutative flux. 

One immediate consequence is that the probe limit may no longer be a valid description of 
the physics. On the macroscopic side, this is due to the addition of U{1) flux through the 5^, 
which effectively contributes a large mass correction to the theory at large n. Therefore, it is 
preferable to consider theories with finite n in order to minimize the effect of back-reaction. This 
is a non-trivial problem in general, with the only conjectured solution being that of the S0{3) 
scalar representation, which has been shown to lead to interesting physics in its own right. 

Given that the finite n prescription and the large n prescription are essentially decoupled 
from one another, we have endeavored in section 2 to show that the finite n action does indeed 
converge to that at large n in a certain limit. This is the first time that such an argument 
has been presented in the literature. How does one then interpret the finite n theory from the 
macroscopic side? Although we have not attempted to answer such a question, we believe that 
the correct interpretation is that of a wrapped D5-brane which now carries non-commutative 
U{1) flux on the world-volume through the introduction of a star product. Thus our D3-brane 
model should be dual to a class of non-commutative field theories which have been explored 
in the literature \60\ I61j . It would be interesting to explore this in more detail since this is a 
different theory to the /t-deformed models and therefore evades the torsional interpretation [6l] 
of the Snyder algebra proposed in [63j, which seems to rule out many non-commutative field 
theory models. 

The above discussion applies for the relativistic theory, and one may be interested in the 
preservation (or not) of non-Abelian physics in the non-relativistic expansion. This was explored 
in sections 3 and 4 of the paper. In section 3, we considered various consequences for the infla- 
tionary scenario when an arbitrary number of coincident branes are present and demonstrated 
explicitly how the slow roll parameters are modified by the non-Abelian structure. We also 
considered a particular small-field inflationary model where the branes are moving out of an 
AdS-type throat. We found that observational bounds on the scalar spectral index can in turn 
constrain the fluxes in the supergravity theory. 

In section 4, we discussed the case of n = 2 in an arbitrary background in both the relativistic 
and non-relativistic limits. In the latter case, we found that the level of non-Gaussianity in the 
primordial curvature perturbation is indeed suppressed, although corrected from the canonical 
slow roll models. We concluded by considering the flow trajectories for such a configuration and 
compared these to those of the standard, single brane DBI scenario. 

We have not addressed the question of reheating in this scenario and this is an important 
topic to consider. Neither have we attempted to identify the location of the standard model. 
On the other hand, we have argued from the field theory perspective that a U{n) field-theory 
with scalars transforming under S0{3)n is dual to a U{1) field-theory with scalars transforming 
under another U{1). This duality appears to be true for all n. It is tempting to suggest that 
one could then identify standard model- type states within a theory (with additional gauge 
fields), after symmetry breaking. This could be a particularly useful description for transferring 
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inflationary energy into the particle sector. 
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